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Abstract 

We present recently obtained results in the theory of pseudoana- 
lytic functions and its applications to elliptic second-order equations. 
The operator (divpgrad-|-(7) with p and q being real valued functions 
is factorized with the aid of Vekua type operators of a special form 
and as a consequence the elliptic equation 

(divpgrad = 0, (1) 

reduces to a homogeneous Vekua equation describing generalized an- 
alytic (or pseudoanalytic) functions. As a tool for solving the Vekua 
equation we use the theory of Taylor and Laurent series in formal 
powers for pseudoanalytic functions developed by L. Bers. The series 
possess many important properties of the usual analytic power se- 
ries. Their applications until recently were limited mainly because of 
the impossibility of their explicit construction in a general situation. 
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We obtain an algorithm which in a really broad range of practical 
applications allows us to construct the formal powers and hence the 
pseudoanalytic Taylor series in explicit form precisely for the Vekua 
equation related to equation (1). In other words, in a bounded domain 
this gives us a complete (in C-norm) system of exact solutions of (1). 



The foundations of pseudoanalytic function theory have been created by a 
considerable number of mathematicians among which Lipman Bers and Ilya 
Vekua played the most prominent role. In the works of I. Vekua and many 
other researchers pseudoanalytic functions are called generalized analytic. 
Nevertheless in the present work we use the term "pseudoanalytic" in order 
to emphasize the fact that we mainly use the part of the theory developed 
by L. Bers and his collaborators. 

In the recent works of the author [29], [30] and [32] a close connection 
between the second-order elliptic equation 



and a Vekua equation of a special form was presented. This connection 
is a direct generalization of a relation which exists between harmonic and 
analytic functions. The special form of the arising Vekua equation (which 
we call the main Vekua equation) allows us to apply well developed methods 
of pseudoanalytic function theory ([4], [8], [17], [41], [43], [44] and others) 
and of p-analytic function theory [39] to the analysis of the corresponding 
second-order equations. In this work we restrict ourselves to the development 
of the theory of series in formal powers for the main Vekua equation and as 
a consequence of the corresponding theory for equation (2). 

Formal powers were defined by L. Bers (see [8]) and represent a gen- 
eralization of the usual powers {z — Zq)" which play a crucial role in the 
one-dimensional complex analysis. As their name reveals formal powers in 
general are not powers. They behave as {z — zq)'^ only locally, near the center, 
and in fact can be complex functions of a quite arbitrary nature. Neverthe- 
less they are solutions of a corresponding Vekua equation and under quite 
general conditions represent a complete system of its solutions in the same 
sense as any analytic function under quite general and well known conditions 
can be represented as its normally convergent Taylor series. 
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Introduction 



(divpgrad -|-g)ti = 



(2) 
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The main result of the present work is a procedure for exphcit construc- 
tion of formal powers corresponding to the main Vckua equation in a very 
general situation. From the relation of the main Vckua equation to equation 
(2) we obtain that under quite general conditions we are able to construct ex- 
phcitly a complete system of solutions of (2). More precisely, let us consider, 
e.g., the conductivity equation 

div (pgradti) = 0. 

Our result then gives us the possibility to construct explicitly a complete 
system of solutions of this equation if p has the form 

p = $((^)vl/(^) (3) 

where ip) is any orthogonal coordinate system, $ and ^ are arbitrary 
positive differentiable functions. 

In the case of the stationary Schrodingcr equation 

(-A + q)u^O (4) 

in order to construct a complete system of solutions explicitly we need a 
particular solution of this equation of the form (3). Note that before this 
result has been obtained the knowledge of one particular solution of a second- 
order equation in two dimensions like (4) had not given much information 
about the general solution. Now one particular solution of (4) is a generator of 
a complete system of solutions of (4) which in a sense and for many purposes 
represents the general solution of the equation. 

The present work is an introduction to this new method and includes 
explanation of the theory behind it and some first examples of application. 
The paper is organized as follows. In section 2 we present the already known 
results from pseudoanalytic function theory concerning mainly the theory of 
series in formal powers. In section 3 we give our recent results on the rela- 
tionship between equation (2) and the main Vekua equation. In section 4 we 
explain how a complete system of solutions of equation (2) is obtained from a 
complete system of pseudoanalytic formal powers. Section 5 is dedicated to a 
classical but unfortunately not widely known relation between orthogonal co- 
ordinate systems in the plane and complex analytic functions. In section 6 we 
first explain and then prove the main new result of this work which opens the 
way for explicit construction of formal powers for the main Vekua equation. 
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In section 7 we explain how this construction gives us complete systems of 
solutions of second-order elliptic equations. We present some examples which 
were calculated using the Matlab tool for symbolic calculation. The choice of 
the examples was inspired by the desire to show another interesting relation 
between our work and the theory of so called i/-regular functions (solutions 
of Vekua equations with constant coefficients) started by R. J. Duffin in [18] 
and [19]. We show that the complete systems of solutions of Yukawa and 
Helmholtz equations obtained in [11] represent an example of a complete 
system of solutions for a second-order elliptic equation constructed on the 
base of pseudoanalytic formal powers. 

Note that some of the results presented in this work are sufficiently simple 
that may serve as a nice and useful addition to a standard complex analysis 
course, for example, the construction of conjugate mctaharmonic functions 
(see subsection 3.3). The author has already done this didactical experiment 
successfully. At the same time it is clear that very close to these simple 
things there are many open and not so simple problems. First of all, the 
theory of pseudoanalytic (or generalized analytic) functions is very rich and 
meanwhile some of the well known results from this theory like, e.g., the 
Cauchy integral theorem and the Morera theorem have been already trans- 
ferred to the corresponding second-order elliptic equation (see [29]), many 
others still wait for their application. For example, it is of great interest 
and importance to construct not only positive formal powers (as we do in 
the present work) but also the negative formal powers for the main Vekua 
equation and hence the Cauchy kernel and integral representations for solu- 
tions of second-order equations. Moreover, in the present work we consider 
equation (2) with real valued coefficients. Nevertheless as was shown in [32] 
the same scheme works in the case of complex valued coefficients in (2), but 
the corresponding main Vekua equation is bicomplex. The theory of bicom- 
plex Vekua equations (which are closely related to the Dirac equation with 
electromagnetic potential) we started to develop in [12], but it is clear that it 
is much more complicated than the theory of complex Vekua equations and 
represents a separate and important challenge. 

The development in three or more dimensions of a theory comparable with 
the theory of pseudoanalytic functions in the plane was an object of study 
of many researchers (see, e.g., [36]). In general, the results obtained in this 
direction are much less complete than their two-dimensional counterparts. In 
[32] a quaternionic generalization of the main Vekua equation was obtained, 
an equation which clearly preserves the basic important properties which 
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made the present work possible. It is interesting to develop the L. Bers 
theory for this main quaternionic Vekua equation which seems to be possible 
due to its quite special form. 

Finally, behind the results presented here there is an essential and in a 
sense universal idea of factorization. The main Vekua equation (as we show 
in section 3) is a result of a factorization of a corresponding second-order 
operator. This kind of factorization is studied not only in the case of elliptic 
operators ([5], [7], [27], [28], [29], [32]) but as well in the case of operators of 
other types ([6] and [13]). Thus, it seems quite possible that similar results 
can be obtained, e.g., not only for a stationary Schrodinger equation (as in 
the present work) but also for a time-dependent Schrodinger equation. 

2 Some definitions and results from pseudo- 
analytic function theory 

This section is based on notions and results presented in [8] and [9]. Let fl 
be a domain in R^. Throughout the whole paper we suppose that Q is a 
simply connected domain. 



2.1 Generating pair, derivative and antiderivative 

Definition 1 A pair of complex functions F and G possessing in Q partial 
derivatives with respect to the real variables x and y is said to he a generating 
pair if it satisfies the inequality 

Im(FG') > inn. (5) 

Denote dj = \ {^-§^ + j and dz = \ — i^^ ■ The following expres- 
sions are known as characteristic coefficients of the pair (F, G) 

^ _ FGz — FzG ^ _ FGz — FzG 
"(■^'^^ ~ FG-FG ' ^^'^^ ~ FG-FG ' 

FG-y, — F-,G FGz — F^G 

JD 



FG-FG ' ^■^'^^ FG-FG ' 
where the subindex 'z or z means the application of dz or dz respectively. 
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Every complex function W defined in a subdomain of admits the unique 
representation W = (f)F + ipG where the functions (f) and ip are real valued. 
Thus, the pair {F,G) generalizes the pair which corresponds to usual 
complex analytic function theory. 

Sometimes it is convenient to associate with the function W the function 
cu — <!) + iip. The correspondence between W and cu is one-to-one. 

The (F, G')-derivative W — '^''"fj^ of a function W exists and has the 
form 

W^cf>,F + ij,G ^W,- A^p,G)W - B^F,G)W (6) 

if and only if 

^^F + V^G^O. (7) 
This last equation can be rewritten also in the following form 

= a(F,G)W + b(F,G)W 

which we call a Vekua equation. Solutions of this equation are called {F, G)- 
pseudoanalytic functions. If W is {F, G)-pseudoanalytic, the associated func- 
tion uj is called (F, G)-pseudoanalytic of second kind. 

Remcirk 2 The functions F and G are (F, G)-pseudoanalytic, and F = G = 
0. 

Definition 3 Let (F, G) and (Fi, Gi) - be two generating pairs in Q. (Fi, Gi) 
is called successor of {F,G) and (F, G) is called predecessor of {Fi,Gi) if 

a(Fi,Gi) = (^{F,G) and ^'(Fi,Gi) = --B(f,g)- 

The importance of this definition becomes obvious from the following 
statement. 

Theorem 4 Let W be an (F, G)-pseudoanalytic function and let (Fi, Gi) be 
a successor of (F, G). Then W is an (Fi, Gij-pseudoanalytic function. 

This theorem shows us that to the difference of analytic functions whose 
derivatives are again analytic, the (F, G) -derivatives of pseudoanalytic func- 
tions are in general solutions of another Vekua equation. 
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Definition 5 Let {F, G) be a generating pair. Its adjoint generating pair 
{F, G)* = {F*, G*) is defined by the formulas 

2F 2G 



FG-FG FG-FG 
The (F, G) -integral is defined as follows 

J Wd^F,G)Z = F{zi) Re J G*Wdz + G{zi) Re J F*Wdz 

where F is a rectifiable curve leading from zq to zi. 

liW = 4>F + ipG is an (F, G)-pseudoanalytic function where and i/j are 
real valued functions then 

Wd^F,G)Z = W{z) - <P{zo)F{z) - ij{zo)G{z), (8) 



zo 



and as F —G — 0, this integral is path-independent and represents the (F, G)- 

antiderivative of W. 

A continuous function w defined in a domain Q is called {F, G')-integrable 
if for every closed curve T situated in a simply connected subdomain of fl, 



r 



wd(^F,G)Z = 



or what is the same 



Re j G*wdz + iReJ F*wdz = 0. (9) 

Theorem 6 An {F,G) -derivative W of an {F,G)-pseudoanalytic function 
W is {F,G)-integrable. 

Theorem 7 Let {F,G) be a predecessor of {Fi,Gi). A continuous function 
is (Fi, Gi)-pseudoanalytic if and only if it is (F, G)-integrable. 

Remark 8 It is easy to see that in the case F = 1, G = i, equality (9) turns 
into the Cauchy integral theorem for analytic functions: J^wdz — 0. 
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2.2 Generating sequences and Taylor series in formal 
powers 

In order to introduce the notion of pseudoanalytic derivatives of arbitrary 
order the following definition is necessary. 

Definition 9 A sequence of generating pairs {(iVi, Gm)}, m — 0, ±1, ±2, . . . 
, is called a generating sequence if {F.m+i,Gjn+i) is a successor of {Fm,Gm)- 
If {Fq, Go) — {F, G), we say that {F, G) is embedded in {{F^, Gm)}- 

Theorem 10 Let {F, G) he a generating pair in Jl. Let he a hounded 
domain, Qi C Jl. Then {F, G) can he emhedded in a generating sequence in 
Oi. 

Definition 11 A generating sequence {{Fjn^Grr^} is said to have period 
^ > if {F„i+^, Grn+n) is equivalent to {Fjn, Gm) that is their characteris- 
tic coefficients coincide. 

Let W be an (F, G) -pseudoanalytic function. Using a generating sequence 
in which {F, G) is embedded we can define the higher derivatives of W by 
the recursion formula 

dz 

A generating sequence defines an infinite sequence of Vekua equations. 
If for a given (original) Vekua equation we know not only a corresponding 
generating pair but the whole generating sequence, that is a couple of exact 
and independent solutions for each of the Vekua equations from the infinite 
sequence of equations corresponding to the original one, we are able to con- 
struct an infinite system of solutions of the original Vekua equation as is 
shown in the next definition. Moreover, as we show in the next subsection 
under quite general conditions this infinite system of solutions is complete. 

Definition 12 The formal power Zm {a, zq] z) with center at zo G coef- 
ficient a and exponent is defined as the linear comhination of the genera- 
tors Fm, Gm with real constant coefficients \, fj, chosen so that \Fm{zQ) + 
^Gmi^o) = a. The formal powers with exponents n — 1,2, . . . are defined hy 
the recursion formula 

Z^^+'\a, zo; z)^{n + 1) f zi"ii(a, ^o; C)^^(f.,g.)C- (10) 

J zo 
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This definition impfies the following properties. 

1. Zm\ci, Zq] z) is an (F^, G^j-pseudoanalytic function of z. 

2. If a' and a" are real constants, then Zm\a!-\-ia" , Zq]z) — a'Zm\l, Zo;z)+ 
a"zln\i,zo;z). 



3. The formal powers satisfy the differential relations 

nZ^^P{a,Zo;z). 



d{Fm,Gm)Zm\a,zo;z) 



dz 

4. The asymptotic formulas 

Z^\a,zo;z) ^ a{z - zo)"", z ^ zq (11) 

hold. 

Assume now that 

oo 

iy(^) = 5^Z(")(a,^o;^) (12) 

n=0 

where the absence of the subindex m means that all the formal powers corre- 
spond to the same generating pair (F, G), and the series converges uniformly 
in some neighborhood of zq. It can be shown that the uniform limit of pseudo- 
analytic functions is pseudoanalytic, and that a uniformly convergent series 
of (F, G)-pseudoanalytic functions can be (F, G')-differentiated term by term. 
Hence the function W in (12) is (F, G) -pseudoanalytic and its rth derivative 
admits the expansion 

oo 

lyM(z) = J] n(n - 1) • • • (n - r + l)Zl^-'\an, Zo] z). 

n=r 

Prom this the Taylor formulas for the coefficients are obtained 

Qn = j . (13) 
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Definition 13 Let W{z) be a given {F,G)-pseudoanalytic function defined 
for small values of \z — Zo\. The series 

oc 

5^Z(«)(a,^o;^) (14) 

n=0 

with the coefficients given by (13) is called the Taylor series of W at Zq, 
formed with formal powers. 

The Taylor series always represents the function asymptotically: 

N 

W{z) - J2 ^^"Ha, ^o; z)=0 (\z - Zo\^^') , z^zo, (15) 

n=0 

for all N. This implies (since a pseudoanalytic function can not have a zero 
of arbitrarily high order without vanishing identically) that the sequence of 
derivatives {^^["^(-Zo)} determines the function W uniquely. 

If the series (14) converges uniformly in a neighborhood of Zq, it converges 
to the function W. 

2.3 Convergence theorems 

The statements given in this subsection were obtained by L. Bers [8], [10] 
and S. Agmon and L. Bers [1]. 

Theorem 14 [8] The formal Taylor expansion (14) of a pseudoanalytic func- 
tion in formal powers defined by a periodic generating sequence converges in 
some neighborhood of the center. 

This theorem means only a local completeness of the system of formal 
powers. The following definition due to L. Bers describes the case when cor- 
responding formal powers represent a globally complete system of solutions 
of a Vekua equation much as in the case of usual powers of the variable z 
and the Cauchy-Riemann equation. 

Definition 15 [8] A generating pair {F, G) is called complete if these func- 
tions are defined and satisfy the Holder condition for all finite values of z, the 
limits F{oo), G{oo) exist, Im(F(oo)G'(oo)) > 0, and the functions F{l/z), 
G{l/z) also satisfy the Holder condition. A complete generating pair is 
called normalized if F{oo) = 1, G{oo) — i. 
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A generating pair equivalent to a complete one is complete, and every 
complete generating pair is equivalent to a uniquely determined normalized 
pair. The adjoint of a complete (normalized) generating pair is complete 
(normalized) . 

Prom now on we assume that {F, G) is a complete normalized generating 
pair. Then much more can be said on the series of corresponding formal pow- 
ers. We limit ourselves to the following completeness results (the expansion 
theorem and Rungc's approximation theorem for pseudoanalytic functions). 

Following [8] we shall say that a sequence of functions Wn converges 
normally in a domain ft if it converges uniformly on every bounded closed 
subdomain of ft. 

Theorem 16 Let W be an {F, G) -pseudoanalytic function defined for \z — zq\ < 
R. Then it admits a unique expansion of the form W{z) = Yl'^=o Z^'"\am zq] z) 
which converges normally for \z — zq\ < 6R, where 9 is a positive constant 
depending on the generating sequence. 

The first version of this theorem was proved in [1]. We follow here [10]. 

Remark 17 Necessary and sufficient conditions for the relation 9 = 1 are, 
unfortunately, not known. However, in [10] the following sufficient conditions 
for the case when the generators {F, G) possess partial derivatives are given. 
One such condition reads: 



for some < £ < 1. 

Theorem 18 [10] A pseudoanalytic function defined in a simply connected 
domain can he expanded into a normally convergent series of formal polyno- 
mials (linear combinations of formal powers with positive exponents). 

Remark 19 This theorem admits a direct generalization onto the case of a 
multiply connected domain (see [10]). 




for some £ > 0. Another condition is 
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In posterior works [24], [37], [20] and others deep results on interpolation 
and on the degree of approximation by pseudopolynomials were obtained. 
For example. 

Theorem 20 [37] Let W he a pseudoanalytic function in a domain bounded 
by a Jordan curve and satisfy the Holder condition on dQ with the exponent 
a (0 < a < 1). Then for any e > and any natural n there exists a 
pseudopolynomial of order n satisfying the inequality 

\W{z)- Pn{z)\<^^^ for any z en 
where the constant does not depend on n, but only on e. 

3 Solutions of second order elliptic equations 
as real components of complex pseudoana- 
lytic functions 

3.1 Factorization of the stationary Schrodinger opera- 
tor 

It is well known that if /o is a nonvanishing particular solution of the one- 
dimensional stationary Schrodinger equation 

then the Schrodinger operator can be factorized as follows 

We start with a generalization of this result onto a two-dimensional situation. 
Consider the equation 

(_A + i.)/ = (16) 

in some domain Q C R^, where A = ^ + i/ and / are real valued 
functions. We assume that / is a twice continuously differentiable function. 
By C we denote the complex conjugation operator. 
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Theorem 21 [30] Let f he a positive in Vt particular solution of (16). Then 
for any real valued function (/? e C^(0) the following equalities hold 



(17) 

Proof. Consider 

% + f C) (9. - ^C) ^ = iA^ - - a. (^) ^ 

= i(A^ - = i (A - ^. (18) 

Thus, we have the first equality in (17). Now application of C to both sides 
of (18) gives us the second equality in (17). ■ 

The operator dz — y /, where / is the identity operator, can be represented 
in the form 

Let us introduce the notation P = fdzf"^!- Due to Theorem 21, if / is a 
positive solution of (16), the operator P transforms real valued solutions of 
(16) into solutions of the Vekua equation 

dz + ^C^w^O. (19) 

Note that the operator dz applied to a real valued function ip can be 
regarded as a kind of gradient, and if we know that = $ in a whole 
complex plane or in a convex domain, where $ = $1 +i$2 is a given complex 
valued function such that its real part $1 and imaginary part $2 satisfy the 
equation 

dy^,+d,^2^0, (20) 

then we can reconstruct (p up to an arbitrary real constant c in the following 
way 

cp{x, y)^2( [ $i(r^, y)drj - f $2(^0, 0^^^) + c (21) 
where (a;o,yo) is an arbitrary fixed point in the domain of interest. 
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By A we denote the integral operator in (21): 



+ C. 



Note that formula (21) can be easily extended to any simply connected do- 
main by considering the integral along an arbitrary rectifiable curve T leading 
from (a;o,2/o) to {x,y) 



(I 



(fi{x, y) = 2 / ^idx - ^2dy + c. 



Thus if $ satisfies (20), there exists a family of real valued functions </? such 
that (9^(^ = $, given by the formula (/? = 

In a similar way we define the operator A corresponding to dj'. 

Z[$](a;,2/) = 2( /" ^i{7],y)d7]+ T $2(0:0, Oc^C) + c 

\Jxo Jyo / 

Consider the operator S — fAf"^! apphcable to any complex valued 
function w such that $ = f~^w satisfies condition (20). Then it is clear that 
for such w we have that PSw — w. 

Proposition 22 [31] Let f be a positive particular solution of (16) and w 
he a solution of (19). Then the real valued function g — Sw is a solution of 
(16). 

Proof. First of all let us check that the function ^ — w/ f satisfies (20). Let 
= Rett; and v — Imw. Consider 

9,$i + 9,$2 = y (^{dyu + d,v) - i^u + ^v)^ . (22) 
Note that equation (19) is equivalent to the system 

O O 9xf dyf ^ ^ dyf dxf 

OxU - OyV = — —u + -j-'"^ ^y"^ + (^xV = + -j-v 

from which we obtain that expression (22) is zero. Thus the function $ 
satisfies (20) and hence the real valued function (p — A[w/ f] is well defined 
and satisfies the equation d^ifi — w/ f. 
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Consider the expression 



d^d,{Sw) = d^i^{dJ)A[j]+w 

= Q A/) A[j] + (dj) ckAij] - ^w. (23) 



For the expression we have 



where the following observation was used 



M + iw / r f u{r],y)\ v{xo,y)\ 



Thus substitution of (24) into (23) gives us the equality 



Proposition 23 [31] Let g he a real valued solution of (16). Then 

SPg^g + cf 
where c is an arbitrary real constant. 
Proof. Consider 

SPg^fAd,[^]^f(^ + c)^g + cf. 

■ 

Theorem 21 together with Proposition 22 show us that equation (16) is 
equivalent to the Vekua equation (19) in the following sense. Every solution 
of one of these equations can be transformed into a solution of the other 
equation and vice versa. 
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3.2 Factorization of the operator divpgrad+g. 

The following statement is known in a form of a substitution (see, e.g., [38]). 
Here we formulate it as an operational relation. 

Proposition 24 Let p and q be complex valued functions, p e C^{fl) and 
p in Q. Then 

divpgrad+g = p^^'^{A — r)p^^^ in (25) 

where 

_ q 

pl/2 p 

Proof. The easily verified relation 

A„l/2 

divpgrad = p'/\A - -jjr)p'^^ (26) 

is well known (see, e.g., [42]). Adding to both sides of (26) the term q (and 
representing it on the right-hand side as p^^^ (q/p)?^^"^) gives us (25). ■ 

Theorem 25 [32] Let p and q be real valued functions, p £ C'^[fl) andp ^ 
in fl, Uq be a positive particular solution of the equation 

{div p grad. +q)u — inQ. (27) 

Then for any real valued continuously twice differentiable function tp the fol- 
lowing equality holds 

^(divpgrad+?)(^ = pV2 (^Q^ + fjc^ _ kc^ ^1/2^^ (28) 

where 

f = p"W (29) 

Proof. This is based on (17). Prom (25) we have that if Uq is a solution of 
(27) then the function (29) is a solution of the equation 

(A-r)/ = 0. (30) 

Then combining (25) and (17) we obtain (28). ■ 
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Remcirk 26 According to (26), A — r — f ^ div grad / ^ where f is a 
solution of (30). Then from (25) we have 

div p grad +g = p^^V"^ div grad (31) 

Taking into account (29) we obtain 

divpgrad+g = Uq^ div pu^ graduQ^ in Q. 

Remark 27 Let q = 0. Then Uq can be chosen as Uq = 1. Hence (28) gives 
us the equality 

idiv(pgrad^) = p'/' (a. + (9. - (p'/V). 

In what follows we suppose that in VL there exists a positive particular 
solution of (27) which we denote by Uq. 

Let / be a real function of x and y. Consider the Vekua equation 

= yF in n. (32) 

This equation plays a crucial role in all that follows hence we will call it the 
main Vekua equation. 

Denote Wi = ReW and W2 = Im W. 

Remcirk 28 [30j Equation (32) can be written as follows 

fd^{f-^W^) + if-^UfW2) = 0. (33) 

Theorem 29 [32] Let W ^ Wi + iW2 be a solution of (32). Then U = 
f~^Wi is a solution of the conductivity equation 

div(/^VC/) = inn, (34) 

and V — fW2 is a solution of the associated conductivity equation 

div(/-VK) = zn fi, (35) 
the function Wi is a solution of the stationary Schrddinger equation 

-AWi + nWi = inn (36) 

with ri = Af / f, andW2 is a solution of the associated stationary Schrddinger 
equation 

-AW2 + r2W2 = inn (37) 
where r2 = 2{Vff/f - r, and {Vff - /J + f^. 
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Proof. To prove the first part of tlie tlieorem we use tlie form of equation 
(32) given in Remark 28. Multiplying (33) by / and applying gives 

from where we have that ^Q{dz{pdz{f~^Wi))) — which is equivalent to 

(34) where U = f'^Wi. 

Multiplying (33) by and applying dz gives 

(f-'W,) + (/-%(/W^2)) = 

from where we have that Re{dz (/~^^ (/W^2))) = which is equivalent to 

(35) where V = fW2. 
From (26) we have 

(A - n) w, = divifv (f-'w,)). 

Hence from the just proven equation (34) we obtain that Wi is a solution of 

(36) . 

In order to obtain equation (37) for W2 it should be noticed that 
/div(/-'V(/W^2)) = (A - rs) W2. 



Remcirk 30 Observe that the pair of functions 

F = f and G = j (38) 

is a generating pair for (32). This allows us to rewrite (32) in the form of 
an equation for pseudoanalytic functions of second kind ( equation (7)) 

Vzf + ipzj = 0, (39) 

where (p and ip are real valued functions. If (p and ip satisfy (39) then W — 
(pf + is a solution of (32) and vice versa. 
Denote w — (f + ipi. Then from (39) we have 

{w + w)zf +{w- w)zj = 0, 
18 



which is equivalent to the equation 

1- P 

The relation between (40) and (34), (35) was observed in [3] and resulted to 
be essential for solving the Calderon problem in the plane. 

Theorem 31 [32] Let W = Wi + iW2 be a solution of (32). Assume that 
f — p^/'^uq, where uq is a positive solution of (27) in Q. Then u — p~^^'^Wi 
is a solution of (27) in Q, and v — p^^^W2 is a solution of the equation 

(div - grad+gi)v = in fl, (41) 
P 

where 

,,..l(l^J^,Y^) + 2(^)']. (42) 

P \P \ P Uo / \ Uo J J 

Proof. According to theorem 29, the function f~^Wi is a solution of (34). 
Prom (31) we have that 

(divpgrad+g) (p-^/^l^i) = /-Miv(/2V {f-^Wi)) 

from which we obtain that u = p~^^'^Wi is a solution of (27). 

In order to obtain the second assertion of the theorem, let us show that 

pi/2(div^grad+gi)(pV2^) = /div(/"2v(/^)) 
for any real valued (p e C^(Q). According to (26), 

/ div(/-2V(/(^)) = (^A - ^) ^=iA- r2) ^. 

Straightforward calculation gives us the following equality 

A/-^ _ 3 fVpY lAp , /Vp Vuo\ Auo , ^ f^uo^ ^ 



/I 4:\pJ 2 p \ P Uo / Uo \ Uo 
From the condition that uq is a solution of (27) we obtain the equality 

Auo _ 5 _|_ /^P 



Uq p \ p Uo 
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Thus, 



/-I 4 \ p y 2 p \ P ' Uq / p \ uo 
Notice that 

Ap-V^ _ 3 /^Vpy _ 1 Ap 



Then 



2 



/I p \ P ' Mo / p V Wo 

Now taking gi in the form (42) we obtain the result from (25). ■ 

3.3 Conjugate metaharmonic functions 

Theorems 29 and 31 show us that as much as real and imaginary parts of 
a complex analytic function are harmonic functions, the real and imaginary 
parts of a solution of the main Vekua equation (32) are solutions of asso- 
ciated stationary Schrodinger equations being also related to conductivity 
equations as well as to more general elliptic equations (27) and (41). The 
following natural question arises then. We know that given an arbitrary real 
valued harmonic function in a simply connected domain, a conjugate har- 
monic function can be constructed explicitly such that the obtained couple 
of harmonic functions represent the real and imaginary parts of a complex 
analytic function. What is the corresponding more general fact for solutions 
of associated stationary Schrodinger equations (which we slightly generaliz- 
ing the definition of I. N. Vekua call metaharmonic functions) and of other 
aforementioned elliptic equations. The precise result for the Schrodinger 
equations is given in the following theorem. 

Theorem 32 [30] Let W\ he a real valued solution of (36) in a simply con- 
nected domain Q. Then the real valued function W2, solution of (37) such 
that W ~ Wi + iW2 is a solution of (32), is constructed according to the 
formula 

W2 = f-'A{^f'^,{f-'W^)). (43) 

Given a solution W2 of (37), the corresponding solution Wi of (36) such 
that W — Wi + iW2 is a solution of (32), is constructed as follows 

Wi = -fA{if-^ch{fW2)). (44) 
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Proof. Consider equation (32). Let W — (j)f + iip / f he its solution. Then 
the equation 

^-tf(h = (45) 

is vaUd. Note that if Wi = Re then (p — M^i//- Given 0, ip is easily found 
from (45): 

It can be verified that the expression A{if'^4>z) makes sense, that is dx{p(t>x) + 

dyif4>y) = 0. 

By theorem 29 the function W2 — f ^t/j is a solution of (37). Thus 
we obtain (43). Let us notice that as the operator A reconstructs the real 
function up to an arbitrary real constant, the function 14^2 in the formula (43) 
is uniquely determined up to an additive term cf~^ where c is an arbitrary 
real constant. 

Equation (44) is proved in a similar way. ■ 

Remark 33 When in (36) ri = and f = 1, equalities (43) and (44) turn 
into the well known formulas in complex analysis for constructing conjugate 
harmonic functions. 

Corollary 34 [32] Let U be a solution of (34). Then a solution V of (35) 
such that 

W = fU + if-^ 
is a solution of (32), is constructed according to the formula 

V^A{if^U^). (46) 

Conversely, given a solution V of (35), the corresponding solution U of (34) 
can he constructed as follows: 

U ^ -A{if-^V^). 

Proof. Consists in substitution of Wi = fU and of W2 — f''^V into (43) 
and (44). ■ 

Corollary 35 [32] Let f = p^^'^uo, where uq is a positive solution of (27) in 
a simply connected domain ft and u be a solution of (27). Then a solution 
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V of (41) with qi defined by (42) such that W = p^^'^u + ip ^^^v is a solution 
of (32), is constructed according to the formula 

V = UQ^A{ipuld^{uQ^u)). 

Let V be a solution of (4-1), then the corresponding solution u of (27) such 
that W — p^^'^u + %p~^l'^v is a solution of (32), is constructed according to 
the formula 

u = —UQA{ip^^UQ'^dz{uQv)). 

Proof. Consists in substitution of / = p^^^Mq, Wi = p^l'^u and = p~^^'^v 
into (43) and (44). ■ 



3.4 The main Vekua equation 

The results of this section show us that the theory of the eUiptic equation 

(divpgrad+g)w = 
is closely related to the equation (32): 

= jW. (47) 

As was pointed out in remark 30 the pair of functions: F = f and G = j 
is a generating pair for this equation. Then the corresponding characteristic 
coefficients ^(f,g) and B(^f^g) have the form 

^(F,G) = 0, B(F,G) = y, 

and the (F, G) -derivative according to (6) is defined as follows 

Comparing B(^p^G) with the coefficient in (19) and due to Theorem 4 we 
obtain the following statement. 

Proposition 36 Let W be a solution of (4V- Then its {F,G)- derivative, 
the function w — W is a solution of (19). 
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This result can be verified also by a direct substitution. 
According to (8) and taking into account that 



F* = -if and G* = l/f, 



the {F, G)-antiderivative has the form 

£ w{Od^F,G)C = f{z) Re £ - Re £ tf{C)w{OdC 



f{z) Re /' ^dC + -f- Im r /(C)«^(C)c^C, (48) 



and we obtain the following statement. 

Proposition 37 Let w be a solution of (19). Then the function 

W{z)= fw{Qd^F,G)C 
J zo 

is a solution of (47). 



3.5 p-analytic functions 

Definition 38 A function ^ = u + iv of a complex variable z = x + iy is 
said to be p- analytic in some domain iff 

Ux = -Vy, Uy — — Vx in Q, (49) 

where p is a given positive function of x and y which is supposed to be con- 
tinuously differentiable. 

The theory of p-analytic functions was presented in [39] . p-analytic func- 
tions in a certain sense represent a subclass of generalized analytic (or pseu- 
doanalytic) functions studied by L. Bcrs [8, 9] and 1. N. Vekua [43] and it 
should be noticed that this subclass preserves some important properties of 
usual analytic functions which are not preserved by a too ample class of 
generalized analytic functions (corresponding details can be found in [39]). 
p-analytic functions found applications in elasticity theory (sec, e.g., [2, 22]) 
and in axisymmetric problems of hydrodynamics (see, e.g., [21], [45]). 
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As was pointed out in remark 30 the function W = ^Z+i-^// is a solution 
of the main Vekua equation (47) if and only if (f) and tjj satisfy the equation 
(39) which is equivalent to the system 

1 1 

In other words is a solution of the main Vekua equation iff its correspond- 
ing pseudoanalytic function of second kind is /^-analytic. Thus, we obtain 
the following connection between the stationary Schrodinger equation and 
the system defining p-analytic functions. 

Theorem 39 Let f be a positive solution of the equation 

-Ag + pg = (50) 

where u is a real valued function and let Wi be another real valued solution 
of this equation. Then the function $ = Wi/ f + ifW2, where W2 is defined 
by (43) is an f^-analytic function, and vice versa, let $ be an f^-analytic 
function then the function Wi = /Re$ is a solution of (50). 

The following relation between solutions of the conductivity equation and 
p-analytic functions is valid also. 

Theorem 40 Let f be a positive continuously differentiable function in a 
domain and let U be a real valued solution of the equation 

div{f^VU) =0 mn. (51) 

Then the function ^ = U + iV is f'^-analytic in Q, where V is defined by 
(46), and vice versa, let $ be f^-analytic in VL then U — Kq<^ is a solution 
of (51). 

Thus, solutions of the stationary Schrodinger equation and of the con- 
ductivity equation can be converted into p-analytic functions and vice versa. 
In some cases this relation leads to a simplification of a part of p-analytic 
function theory. 

Example 41 [31] A considerable part of bibliography dedicated to p-analytic 
functions consists of studying the case p — , where A; e R (see, e.g., [16], 
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[25], [39]). Let us see what is the form of the corresponding Schrodinger 
equation. For this we should calculate the potential u in (50) when f — x^l"^ . 
It is easy to see that 

The Schrodinger equation with this potential is well studied. Separation of 
variables leads us to the equation 

X"{x) + [^-^ - X{x) = 0, (53) 

where 0^ is the separation constant and a — {k — l)/2. The function 

X{x) = \/xZa{l3x) 

is a solution of (53) (see [23, 8.491]) where denotes any cylindric function 
of order a (Bess el functions of first or second kind). Thus the study of x^- 
analytic functions reduces to the Schrodinger equation (16) with u defined by 
(52) which in its turn after having separated variables reduces to a kind of 
Bessel equation (53). 

Example 42 In the work [26] boundary value problems for p- analytic func- 
tions with p — xj [x?' + y^) were studied. Considering 



x^ + 

we see that this function is a solution of the Schrodinger equation (16) with 
V having the form 

_ 1 

that is we obtain again the potential of the form (52) where k = 1 and as was 
shown in the previous Remark the study of corresponding p-analytic functions 
in a sense reduces to the Bessel equation (53). 



4 Complete systems of solutions for second- 
order equations 

In what follows let us suppose that the real valued function / is defined in a 
somewhat bigger domain with a sufficiently smooth boundary. Then we 
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change the function / for z G and continue it over the whole plane in 

such a way that / = 1 for large \z\ (sec [10]). In this way the generating pair 
{F,G) = (f,i/f) becomes complete and normalized. 

Then the following statements are direct corollaries of the relations estab- 
lished in section 3 between pseudoanalytic functions (solutions of (32)) and 
solutions of second-order elliptic equations, and of the convergence theorems 
from subsection 2.3. 

Definition 43 Let u{z) be a given solution of the equation (27) defined for 
small values of \z — Zo\, and let W{z) be a solution of (32) constructed ac- 
cording corollary 35, such that ReW — p^^'^u. The series 

CXD 

p-V2(^)^ReZ(")(a„,^o;^) 

n=0 

with the coefficients given by (13) is called the Taylor series ofu at zq, formed 
with formal powers. 

Theorem 44 Let u{z) be a solution of (27) defined for \z — Zo\ < R. Then 
it admits a unique expansion of the form 

oo 

u{z) ^p-^/\z)J2^^Z^''\an,zo;z) 

n=0 

which converges normally for \z — Zo\ < R. 

Proof. This is a direct consequence of theorem 16 and remark 17. Both 
necessary conditions in remark 17 are fulfilled for the generating pair (38). 
■ 

Theorem 45 An arbitrary solution of (27) defined in a simply connected 
domain where there exists a positive particular solution Uq can be expanded 
into a normally convergent series of formal polynomials multiplied byp"^/'^. 

Proof. This is a direct corollary of theorem 18. ■ 

More precisely the last theorem has the following meaning. Due to Prop- 
erty 2 of formal powers we have that Z^'^\a, Zq; z) for any Taylor coefficient 
a can be expressed through Z("^(l, zq] z) and Z^^\i, zq; z). Then due to the- 
orem 18 any solution W of (32) can be expanded into a normally convergent 
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series of linear combinations of Z^'^\l, z^; z) and Z^^^i, Zq; z). Consequently, 
any solution of (27) can be expanded into a normally convergent series of 
linear combinations of real parts of ^o; and Z^'^\i, zq; z) multiplied 

by 

Obviously, for solutions of (27) the results on the interpolation and on 
the degree of approximation like, e.g., theorem 20 are also valid. 

Let us stress that theorem 45 gives us the following result. The functions 

{p-''\z) ReZH(l, zo; z), p-''\z) ReZW(i, z,-, z)]°^^^ (54) 

represent a complete system of solutions of (27) in the sense that any solu- 
tion of (27) can be represented by a normally convergent series formed by 
functions (54) in any simply connected domain Q where a positive solution of 
(27) exists. Moreover, as we show in section 6, in many practically interesting 
situations these functions can be constructed explicitly. 

5 A remark on orthogonal coordinate sys- 
tems in a plane 

Orthogonal coordinate systems in a plane are obtained (see [35]) from Carte- 
sian coordinates x, y by means of the relation 

u + iv = $(a; + iy) 

where $ is an arbitrary analytic function. Quite often a transition to more 
general coordinates is useful 

^ and 1] preserve the property of orthogonality. Some examples taken from 
[35] illustrate the point. 

Example 46 Polar coordinates 

u + iv — \n{x + iy), 

u — In V — arctan — . (55) 

X 

Usually the following new coordinates are introduced 

r — — y/x^~-\-y^, ip — v — arctan — . 
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Example 47 Parabolic coordinates 

u + iv 



V2 



u = yr + x, V = \Jr — X. 
More frequently the parabolic coordinates are introduced as follows 

t- 1 2 

Example 48 Elliptic coordinates 

X + iy 



u + tv — arcsm ■ 



a 



Si — S2 , Si + S2 

smu — , cosni> = — 

2a 2a 

where si = \/ {x + a^ + y"^, S2 = ^/ {x — a)^ + y'^. The substitution 

^ — sin u, T] — cosh v 

is frequently used. 

Example 49 Bipolar coordinates 

a + X + iy 

u + iv — in , 

a — X — ly 

, 2ax 2ay 

tsiiahu=— -, temv = — r -. 

a^ + x"^ + a^ — x^ — y^ 

The following substitution is frequently used 

6 Explicit construction of a generating se- 
quence 

We suppose that the function / in the main Vekua equation (32) has the 
following form 

/ = U{u)V{v) (56) 
28 



where u and v represent an orthogonal coordinate system and according to 
the explained in the previous section we assume that $ = u + iv is an analytic 
function of the variable z = x + iy. U and V are arbitrary differentiable 
nonvanishing real valued functions. 

The present section is structured in the following way. First we explain 
how one can naturally arrive at the form of generating sequences in the main 
result of the section, theorem 50 and then we give a rigorous proof of this 
result. 

The first step in the construction of a generating sequence for the main 
Vekua equation (32) is the construction of a generating pair for the equation 
(19) which as was shown in subsection 3.4 is a successor of the main Vekua 
equation. For this one of the possibilities consists in constructing another 
pair of solutions of (32). Then their (F, G)-derivatives will give us solutions 
of (19). Let us consider equation (32) in its equivalent form (39): 

fz + j^^z^O (57) 
and look for a solution in the form 

Then we have the following equation 

(p'{u)u^+ —iIj'{v)v^ = 0. 

Taking into account that $ = u+iv is an analytic function, that is Uz+iVz = 
we observe that (57) is fulfilled ilip'{u) = ip^v)/ f^. Now using (56) we obtain 
that U^{u)V^{v) = ip'{v)/^'{u) and hence 

and V'(^) = J V''{v)dv. 

The corresponding solution Wi of (32) has the form 

Wi^ f -^U{u)V{v) + / V^{v)dv- ' 



U^{u) ' ' ' ' 7 ' ' U{u)V{v) 

Its (F, G)-derivative is obtained according to (6) as follows 

V V V 
Wi = —Uz + i—Vz = —^z- 
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By analogy, we can look for a solution of (57) in the form 

(p = ^{v), il) = ip{u). 
Then we have the following equation 

ip'{v)v^+ —ip'{u)u^ = 0. 

This equation is fulfilled if (p'{v) = —1/V'^{v) and ip'{u) = U^{u). Conse- 
quently, 

^{v) = - J and ^(m) = J U^{u)du. 

The corresponding solution W2 of (32) has the form 



= - y" ^^U{u)V{v) + J U\u)d 



u 



U{u)V{v) 



with its corresponding (F, G)-derivative 

u -U .u ^ 
^ V V V 

Denote Fi = Wi and Gi = W2- It is easy to see that the pair Fi, Gi fulfils 
(5) and by construction satisfies equation (19). Thus, (Fi, Gi) is a successor 
of (F,G) defined by (38). 

The following step is to construct the generating pair {F2, 6*2) . For this we 
should find another pair of solutions of (19) and apply the {Fi, Gi)-derivative 
to them. Consider equation (19) written in the form 

VzFi + i^zGi = 
which in our case can be represented as follows 

V U 

Vzjj + ^ziy^O- (58) 
Again, let us look for a solution in the form 
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Then equation (58) is satisfied if (p'{u) = (y^j^ i})'{v) from where we obtain 

^{u) = J U^{u)du and ijj{v) = J V^{v)dv. 
Then the corresponding solution of (19) has the form 



wi= / U^{u)du^,--^+ / V'{v)dv^ 



'U{u) J ^ ' " V{v) 
Its (Fi, Gi)-derivative is obtained as follows 

= ^^UVu, + i^^UVv^ = UV {^,f . 
Analogously, looking for a solution of (58) in the form 

we obtain that 

dv , , . X f du 



^i^)-- j and ij{u)^ J 



The corresponding solution of (19) has the form 

f dv ^ V(v) f du ^ iUiu) 



V^{v) 'U{u) J U^{u) ' V{v) ■ 
Its (Fi, Gi) -derivative is obtained as follows 

$2^^ ^Ju^ i 2 

We have that the pair F2 — Wi and G2 — W2is a successor of {Fi,Gi). 
Observe that 

(i^i,Gi) = U'^^G 
and 

(F2,G2) = (($.)' F, i^.fG). 

Prom these formulae it is easy to guess a general form of the corresponding 
generating sequence which we give in the following statement. 
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Theorem 50 Let F — U{u)V{v) and G = ^(^jVct;) ^^^'^^ ^ ^'^^ ^ ^'^^ 

arbitrary dijferentiable nonvanishing real valued functions, $ = -u, + iv is an 
analytic function of the variable z = x + iy in Q such that $2 is bounded 
and has no zeros in fl. Then the generating pair {F, G) is embedded in the 
generating sequence {Fm, Gm), m — 0, ±1, ±2, . . .m Q defined as follows 



F^ = {^z)"^ F and Gm = {^z)"^ G for even m 

and 

-F and Gm = ($^)"' U^G for odd m. 



™ 1/2 

Proof. First of all let us show that {Fm, Gm) is a generating pair for m — 
±1, ±2, .... Indeed we have 

lm(FmGm)^lm{\^zf"'FG)>0. 

Consider a(^F^^Gm)- For both m being even or odd we obtain 

We should verify the equality 

hFm,Gm) = -B(Fr^_,,Gm-l) (59) 

Consider first the case of an odd m. A direct calculation gives us 



and 



(Fm-l,Gm-l) - } 



Thus equality (59) is true iff = {b(F,G) —'^U'zjj) is equal to —Bi^f,g)- It is 
easy to see that 

'u' y 



B(F,G) = Uz[ — -i^ 



Consider 



A o U'\ Uz + ivz fU' V ^ U' 

= \f,g) - 2m^— = — —Uz + —Vz - 2uz — 

^z \ U J Uz — IVz \ U V U 
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Thus, equality (59) is proved in the case of m being odd. 
Now let m be even. Then 



and 



B(p^_, n^_A = { = \ B(p.a) - 2—u 



{Fr^-i,Gm-i) = I f= y y^(F,G) - ^Jj^Z 



Equality (59) is valid iff the expression ^b(^F,G) is equal to (^—B(^pG) + '^jf'^z) ■ 
We have 

and from the other side 

U' fU' .V'\ U' fU' .V'\ 

-B^F,G) + 2— li. = -u, 2—u, ^u,\^— + i—j. 

Thus, equality (59) is proved in all cases and the sequence {Fm,Gm), n^- = 
0, ±1, ±2, . . . satisfies the conditions of Definition 9. Therefore it is a gener- 
ating sequence. ■ 

Remcirk 51 This result obviously generalizes the explicit construction of a 
generating sequence in the case when u = x and v — y presented in a number 
of works by L. Bers (see, e.g., [8] and [17]). 

The last theorem opens the way for explicit construction of formal powers 

corresponding to the main Vekua equation (32) in the case when / has the 
form (56) and hence for explicit construction of complete systems of solutions 
for corresponding second-order elliptic equations. We give more details on 
this as well as some examples in the next section. 
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7 Explicit construction of complete systems 
of solutions of second-order elliptic equa- 
tions 

In section 4 we explained how complete systems of solutions of second-order 
elliptic equations can be constructed from systems of formal powers for a 
corresponding main Vekua equation (32). In the preceding section we estab- 
hshed a result which allows us to make this construction possible in the case 
when / in (32) has the form (56). Then formal powers are constructed simply 
by definition 12. The meaning of this result for the stationary Schrodinger 
equation and for the conductivity equation we discuss separately in the next 
two subsections. 



7.1 Explicit construction of complete systems of solu- 
tions for a stationary Schrodinger equation 

Consider the equation 

-A^ + ug^O in Q. (60) 

where v is a. real valued function. In order to start the procedure of construc- 
tion of a complete system of solutions of this equation we need a positive 
particular solution / having the form (56). 

Example 52 Let u in (60) depend on one Cartesian variable: v = v{x). 
Suppose we are given a particular solution f — f{x) of the ordinary differ- 
ential equation 

+ (61) 

It is sufficient for the application of our result from the preceding section for 

constructing the corresponding generating sequence and hence the system of 
formal powers for the main Vekua equation which in this particular case has 
the form 

= —W. 
2/ 

Example 53 A number of works (see, e.g., [I4] and [15]) are dedicated to 
construction (in our terms) of a particular solution for the Schrodinger equa- 
tion with a radially symmetric potential. This solution f is precisely the only 
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necessary ingredient in order to obtain a complete system of formal powers 
and hence of solutions of the Schrodinger equation. Here an important re- 
striction is that the analytic function $2 corresponding to the polar coordinate 
system has the form 1/z and hence our procedure works in any domain not 
including the origin where f is positive. 

Example 54 Consider the Yukawa equation 

(-A + c^)u = (62) 

with c being a real constant. Take the following particular solution of (62): 
f — e'^ . Let us construct the first few corresponding formal powers with 
center at the origin. We have 

(1, 0; z) = e'^y, (i, 0; z) = ie'^^ 



c 



sinh ( q/ ) 2ixsmh.{cy) 



C2 



Z(^) i^, 0; .) = J^^^^ + ^ ( + !) - ' ■ ■ ■ ■ 

It is a simple exercise to verify that indeed the asymptotic formulas (11) hold. 
Now taking real parts of the formal powers we obtain a complete system of 
solutions of the Yukawa equation: 

( \ cy I \ cy /X Sinh(c7/) 

UQ[x,y)^e\ ui{x,y)^xe% U2[x,y) ^ , 

c 

/ N / 2 y\ cy . sinhjcy) 2xsmh{cy) 

U3{x,y)^[x e*'H , U4{x,y) ^ 

\ c/ c c 

Formal powers of higher order can be constructed explicitly using a computer 

system of symbolic calculation. For this particular example (together with 

Maria Rosalia Tenorio) Matlab 6.5 allowed us to obtain analytic expressions 

for the formal powers up to the order ten, that gave us the first twenty one 

functions Uq, . . . , U2o- We used them for a numerical solution of the Dirichlet 

problem for the Yukawa equation with very satisfactory results. For example, 

in the case when fl is a unit disk with centre at the origin, c—1 and u on the 

boundary is equal to ( this test exact solution gave us the worst precision 
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because of its obvious "disparateness" from functions Uq, Ui . . .) the maximal 
error maXzfzn \u{z) — u{z)\ where u is the exact solution and u = Xln=o ^n^n, 
the real constants a„ being found by the collocation method, was of order 
10~^. A very fast convergence of the method was observed. 

Although the numerical method based on the usage of explicitly or nu- 
merically constructed pseudoanalytic formal powers still needs a much more 
detailed analysis these first results show us that it is quite possible that in due 
time and with a further development of symbolic calculation systems it can 
rank high among other numerical approaches. The use of complete systems 
of exact solutions based on pseudoanalytic formal powers has an important 
advantage of universality. The system does not depend on the choice of a 
domain whenever the function f keeps to have no zeros. 

In [11] the following system of solutions of the Yukawa equation (62) was 
obtained from completely different reasonings, 

uq — cosh(cy), ui — xcosh(cy), ^2 = x^cosh(cy) sinh(cy), 

c 

Us — cosh(q/) sinh(cy), 

The same system is easily obtained ( our Matlab program gives symbolic repre- 
sentations of these functions up to the subindex 21) if instead of the particular 
solution f — e"^^ considered in the above example we choose f — cosh(cy). 
In the case of the Helmholtz equation 

(A + c^)u = 

L. R. Bragg and J. W. Dettman constructed the following system 

y 

Uq = cos(cy), Ui = X cos(cy), U2 = x^ cos(cy) sin(cy), 

c 

U3 = x'^ cos(cy) sm(c|/), .... 

This system is also easily obtained as explained above choosing f = cos(cy). 

Thus, the construction of the systems of solutions for the Yukawa and for 
the Helmholtz equation as well as the proof of their completeness in [11] are 
corollaries of a more general theory presented here. 
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In fact, having a positive particular solution / of the form (56) for equa- 
tion (60) allows us to construct explicitly a complete system of solutions of 
(60) in any simply connected bounded domain where the corresponding to a 
chosen coordinate system analytic function $^ is bounded and different from 
zero. 

7.2 Explicit construction of complete systems of solu- 
tions for a conductivity equation 

Consider the equation 

div(/^VC/) = 0. 

Usually in practical applications the function / is positive and equals 1 out- 
side a certain disk. We suppose additionally that / is continuously differen- 
tiable on the whole plane. Under these conditions F — f and G = i/f is 
a complete normalized generating pair and all the convergence and Runge 
type results are valid. Using results of section 6 one can construct explicitly 
a complete system of solutions of the conductivity equation when the func- 
tion / has the form (56) in any domain where the analytic function from 
theorem 50 is bounded and has no zeros. 
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